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On ■ For nonlinear wave equations with a potential term we prove pointwise space- 

. time decay estimates and develop a perturbation theory for small initial data. We 

CN , show that the perturbation series has a positive convergence radius by a method 

which reduces the wave equation to an algebraic one. We demonstrate that already 
first and second perturbation orders, satisfying linear equations, can provide precise 
information about the decay of the full solution to the nonlinear wave equation. In 
a forthcoming publication (part II) we address the issue of optimal decay estimates 
^ , and precise asymptotics under spherical symmetry where the perturbation equations 

p ^! can be solved almost exactly. 

■ I. INTRODUCTION 



O 



It is a well known fact that the presence of a long-range potential term (power-law decay 
■ at spatial infinity) in the wave equation violates the Huygens principle and gives rise to a 
late-time tail in the solution with a power-law decay (both powers are related) flg. It is 
OO . not so well known that nonlinear terms like cause the same effect. We study equations 
^ ' where both these effects are present and give pointwise decay estimates on the solutions. 

Further, we develop perturbation theory for these equations and by its means argue that 
presented estimates give optimal decay rates at late times. A rigorous proof of this fact will 
! appear in a following publication jl] (part II). 

We consider linear and nonlinear wave equations with a potential term of the general 
form 

><| . Uu + Vu = F{u) (1) 

d \ in 3 spatial dimensions, i.e. u : (t, x) E M+ x M'^ = M}^'^ — > M, and solve the initial value 
problem with 

uiO,x) = f{x), dtu{0,x)=g{x). (2) 

First, we construct an iteration scheme and show its convergence in a weighted space-time 
L°°-norm what reproduces the decay estimate from [2| 

C 

\u(t, x) I < --— , , , V(t, x) e 

' ^ ^' - {l + t+\x\){l + \t-\x\\)i-^ ^ ' ^ + 

with q := min(m — l,k,p — 1) provided the potential V and the initial data f,g satisfy 
pointwise bounds 

\fi^)\<T, rV^' |V/(a:)| < „ , \g{x)\< m>3. 

^' ~ (1 + Ixl)"^-! ' ^ ^' ~ (1 |a;|)™ ^' ~ (1 -f- l^l)"^ 
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with small Vq, /o, /i, /2 and the nonlinearity is analytic and satisfies for p > 1 + \f2 

■ul, \v\ < 1. 

Next, we construct a perturbation series representing the solution u and prove its conver- 
gence (with finite convergence radius) in the same weighted space-time L°°-norm. It implies 
pointwise convergence in M^^^ what allows us to control the decay at every perturbation 
order and obtain estimate on the remainder of the perturbation series for any order. Fi- 
nally, if we can show that at some perturbation level our decay estimate is optimal, i.e. we 
know the true asymptotics for late times (what is not very difficult because the perturbation 
equations are linear) then we immediately know the asymptotics of u. It is the same as that 
of the given perturbation order because all higher terms in the perturbation series, summed 
up, are too small to be able to modify the asymptotics. The issue of optimal decay estimates 
and precise asyniptotics compared with numerical results will be addressed in a forthcom- 
ing publication [3] which will be focused on spherical symmetry where the perturbation 
equations can be solved almost exactly. 

The proof of convergence of the perturbation series is essential for justifying the pertur- 
bation scheme as a rigorous approximation and being able to provide exact decay rates. We 
show it by relating the (inverted) wave equation 

u = n-^F{u) - n-\Vu) + £/(/, g) 

(where el{f,g) stands for initial data contribution to the solution of the free wave equation 
□m = 0) to an algebraic equation of a similar form 

W = CF{W) +6W + eD, 

{F is obtained from F by transformation of its Taylor series), which arises from comparison 
of the perturbation schemes for both problems. We make an interesting observation that 
the nonlinear wave equation has a solution u{e) analytic in e, and hence representable by a 
convergent series in e, if the same holds for the solution W{e) of the corresponding algebraic 
equation. The latter, however, is always true when F{u) is analytic at m = what we 
assume. 

Regarding regularity, we can go a safe way and consider only the classical solutions, i.e. 
assume {f,g) G C\R^) x C\R^), V G C\R^) and F G C\R) and obtain u G C2(M^+3)_ 
However, all results remain true also for weak solutions where {f,g) G C^(M^) x C°(M'^), 
V G C°(M^) and F G C°(M) and we have u G C°(M++^), because the lemmas [H which 
constitute the main "engine" of all estimates, preserve the continuity (see ^ for a detailed 
discussion of the weak solutions). 

This paper is organized as follows. It has three main sections addressing the linear wave 
equation with potential, nonlinear wave equation without and with potential, respectively. 
The idea is to develop tools for the simplest, linear problem and then to generalize them 
to the nonlinear situation. Every section has subsections presenting an iterative and a 
perturbative approach to the construction of solutions and a discussion of the optimal decay 
rates. Appendix collects some lemmas used in the proofs, cited from other works. 

Notation 

With the symbol (x) := 1 + |a;| we define spatial and space-time weighted-L°° norms 

Wfh- := ||(a;r/(x)|U.(M3) 
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IklU- := \\{t+ IxD^t-lxiy ^M(t,x)||^oo(Ri_+3) 
of which we will most frequently use 

hhr:^ ■= IK^+ |a;|)^"^M(t,x)||^oc(K^+3). 

Its finiteness guarantees the decay of u like 1/t on the lightcone t ~ |x| and like for 
fixed X as well as for fixed t. Note that functions with compact support in M.^ belong 

to all spaces with any m > 0. 

We introduce the following notation for solutions of the wave equations. Let Jy be a 
linear map from the space of initial data to the space of solutions of the wave equation 
([TJ-Q with F{u) = 0, so that u = Iv{f, g)- For wave equations with a source term and null 
initial data 

Du + Vu = F, u{0,x) = 0, dtu{0,x)=0, 

let's denote the solutions by m = Lv{F), where Ly is a linear map from the space of source 
functions to the space of solutions to the above problem. Note that, due to linearity, the 
solution M of a wave equation with source F and non- vanishing initial data f,g is a. sum of 
these two contributions 

u = LviF) + Iv{f,g). 
Observe that if we put the potential term on the r.h.s. we obtain 



Uu = -Vu + F 



which, treated as a wave equation without potential (on the l.h.s.), is formally solved by 

u = -Lo{Vu) + Lo{F) + Io{f,g). 

Here the solution u appears on both sides what seems to make the formula useless, but it 
will allow us to formulate various iteration schemes, e.g. 

Un+l = -Lo{VUn) + Lo{F{Un)) + /o(/, g) 

for which we will prove convergence in suitable norms. 

Finally, we define constants which arise from estimates proved in [2|, improved in 



Cm '■= niax ( — -, 5 



9 

2(m~ 21 



Cr,,a :— 2 H - + 



■ p-1 ■ q-1 



The latter will be referred to as a bound on the allowed strength Vq of the potential. Our 
purpose is to emphasize that this bound, although not optimal, is finite and not arbitrarily 
small what is crucial when a potential with a given value Vq is studied (like e.g. in the 
Regge- Wheeler equation describing waves on Schwarzschild geometry). 
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II. LINEAR CASE WITH POTENTIAL 

First, we consider a linear wave equation 

Uu + \V{x)u = (3) 

where A > is a small parameter, bounded by some finite constant Cy > (which will be 
defined later). We first show that a standard iteration scheme converges for all A < Cy to 
a solution in Lf^p, i.e. there exists a constant C such that 

C 

\u{t, x) I < —-— V(t, x) G M+ X R3 

with some p > 2 provided the potential V and the initial data /, V/, g are (at least) contin- 
uous and satisfy pointwise bounds 

in^)l<^' ^>2 (4) 

and 

1/(^)1 <(^' \^f(-^\^J^^ 

Then, we show that a perturbation scheme based on expansion in powers of A is, due to 
linearity, equivalent to the iteration scheme and the perturbation series has convergence 
radius Cy- 

As next, we show that the lowest order uq has, in general, a different decay estimate 
than all higher orders, starting from Ui. Finally, we prove that either mq or Ui gives precise 
information about the decay rate of the full solution u. 



A. Iteration 

We define an iteration by 

M_i := 

Un ■= loif, g) - ALo(l^M„_i), = 0, 1, 2, ... 

Then we have the following 

Theorem 1. With /, g and V as above for any m > 3 and k > 2 the sequence Un converges 
(in norm) in L'^p for p = min(A;,m — 1) provided A < Cpl- The limit u := lim„^ooW„ 
satisfies 

with some positive constant C depending only on /o, fi,go, A and k,m. 

This theorem was proved first for classical solutions in [1] and later generalized to weak 
solutions in [2] and stated in a more detailed form, which will be important here. We cite 
the essential part of the proof because some of the presented estimates will be used later. 
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Proof. For g,V f G L'^^ and / G -^^^-i with m > 3, from lemma [T], we get Uq = Io{f,g) G 
-^i^m-i- Next, observe that if m„ G with some p > 1 then 

||(x)''\^m„||lj°p < ||(a;)''V"||L=o||M„||Loc^ = ||m„||l^^ < oo 



and from lemma [2] with F = Vun we get LoiVun) G Lfp when p < k. Because Lf^p-^ C Lf^p^ 
when pi > we get m„+i G with p < min(m — l,k). By induction we obtain u„ G 
for every n = 0, 1, 2, ... with the optimal value p := min(m — l,k). Then, we have 

ll^n+l - ^nlU- = \\\Lo{-V{Un - M„_i))||ljo^ < \Cp^k\\ {x)^V {Un - Mn-l)|U- 

again making use of lemma [2] with {x)^F = —{x)^V{un — Wn-i) ^ -^i^p- For 5 := XCp^k < 1 
the iteration is a contraction in the normed space L'^p. A simple argument shows that the 
sequence u„ is Cauchy. We have 

WUn+l-Unhf;^ < - n_i || Ljo^ = 5"'^'^ \\Ioif , g)\\ L^^ (6) 

and for n' > n 

n'~n—l n'—n—l 
Ikn' -MnlU- < hn+j+l-Un+j\\Lf;^< ^ 5^+"+^ || Jq (/, ^) || 

i=o j=o ' (7) 



< 



1-5 



i|/o(/,^7)||Lj^,. 



This expression can be made arbitrarily small (smaller than any e > 0) for all n,m > M(e). 
Hence, m„ is a Cauchy sequence in Lf^p which is Banach and has a limit u G Lf^p satisfying 

u = Io{f,g)-XLo{Vu). (8) 

This equation is equivalent to the wave equation ([3]) with the initial data ([5]). Finally, we 
find the Lf^^-norm of u 

< \\Ioif, g)\\Lf:^ + X\\Lo{Vu)\\lo.^ < Cmifo + fi + 9o) + XCp^klMif;^, 

thus 

□ 



B. Perturbation series 

Now, we define a perturbation series by 

oo 

n=0 
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and insert into the wave eq. ([3]). It leads to the following perturbation scheme 

□^0 = 0, {vo,Vom = {f,9) ^ Vo = Io{f,9) (9) 

□t;„+i = -■l/t;„, (f„+i,w„+i)(0) = (0,0) Vn+i = -Lo{Vvn) (10) 

Due to linearity of ([3]) it turns out that the partial sums 



k=0 



give the elements u„ obtained above by the iteration technique, so both methods (if they 
work) are equivalent. Theorem [1] implies convergence in Lf^^ for A < C~l with k > 2, m > 3 
and p := mm{k, m — 1). 

From in the proof of Theorem [T] it follows that 

= Yn ^ {Cp,kyVo{f,9)\\L-, 

hence f„ G L'^p for all n > 0. Observe, however, that in the case when m — 1 > k = p 
we have at the lowest order vq = uq a. better decay estimate, namely vq G Lf^_-^ (see 
first line of the proof). The reason that vq decays faster is that its decay comes only from 
initial data and is not influenced by the potential. At all higher orders, Vn{n = 1,2,...) 
contain the contribution from the scattering on the potential and are only in Lf^f.- Since 
u G L'^p = L'^ky expect that all m„ starting from Ui G L'^^ predict qualitatively correct 
asymptotic behaviour of u while the lowest order uq G L'^m-i ^^i^^ ^^i^- This becomes 
especially evident for initial data with compact support, for which uq G Lf^g with arbitrarily 
big g, but Ml, M2, ... G L'^f, 3 u. 

Knowing that the perturbation series converges for some A we can estimate the error 
of the n-th perturbation's order relative to the exact solution by estimating the sum of all 
higher order terms. For the convergent sequence Un we use the relation ([7]) which holds also 
in the limit n' — >■ cxo, Un' — »■ u and gives 

llw-MnlU- < Y3^II^o(/,C/)||l- • 

It provides a pointwise bound on the error 

X) - .„(t, X) I < . ^,^|^|^^,_|^|^,„i V(t, X) G . (11) 



C. Optimal decay estimate 

In this section we sketch a proof how, under some conditions, the optimal decay estimate 
and precise asymptotic behaviour of the solution u can be deduced from the behaviour of 
low order perturbations. This will be studied in more detail in a forthcoming publication 
[sj dealing with spherical symmetry where the lowest perturbation orders can be calculated 
almost explicitly. 

Consider first the case m — 1 > = p, i.e. when the rate of decay of u is dominated by 
scattering on the potential (and not by decay of the initial data). We have = f o ^ 
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and Vn G for n > 1. Below, we show that if the asymptotic behaviour of Vi is such as 
provided by its estimate (i.e. p in the norm is optimal) then theorem [1] gives an optimal 
estimate for u G L'^p with the same decay rate p. Here, we consider only the asymptotics in 
direction of timelike infinity (the case of spatial infinity can be treated similarly). Assume, 
we are able to show (by some explicit calculation, like in [sl) that vi(t,x) = Lq{Vvq) = 
Ci{x)t~'P 7^ for t ^ 1, where Ci(x) is independent on A. The approximation sign means 
that for every small 1] > and every x G there is a TQ{x,r]) > such that for all 
t > To(a;, 77) the relative error is small, i.e. 



vi{t,x) - 



From flTTl) with Ui = Vq + Xvi we have 



ci{x) 



tP 



< T] 



\ci{x)\ 

tp 



(12) 



\uit,x)-vo{t,x) - Xv^{t,x)\ < • =■ Ai(t,x) 



l-Cp,fcA {t+\x\){t-\x\y 



for all (t,x) G . A simple inequality^ 



< 



implies 



;i-C)'^ - i-ac 
1 



1 + 



1-aC 



< 2, VC < l/(2a), a > 1 



(13) 



(^-1^1)^ (1 + t) 



Ml 



\x\ 
T+i 



< 



[i + tY 



for ( := |x|/(l + t) < l/(2g), hence is true for alH > 2g|x|. The error term can be estimated 



Ai(t,x) <2(Cp,fcA)^ 



.2C^-(/o + /i + 5o) 



C- 



A^ 



(l + t)p (l + t)P 

where we have twice used (fT3|l for t > 2{p — l)\x\ and A < 1/(2 Cp^k)- Further, 



~ A^ 

^ < r/A- 



ci(x)| 



;i + t)P ~ ' tP 

provided A is small enough such that A < A^(x) := rjCi{x)/C. Again from ( |T3|) we get 

Co ^ 2co 



\VQ{t,x)\ < 



{t+\x\){t-\x\) 



m-2 



< 



1 + t)™-l 



for all t > 2(m — 2)|x|. Then, 



\vo{t,x)\ < 



2Cn 



CliX 



[1 + 1)™-^ ~ ' tp 

provided t > Ti{x,ri, A) is big enough, such that t"^-^-P > 2co/{ri X |ci(x)|). 



It follows immediately from the Bernoulli's inequality. 
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Finally, we arrive at the statement that for every small rj > and every x G 
M.^, for sufficiently small A < min[A^(a:), 1/(2 Cp^)], and for sufficiently big t > 
max[To(x, r^), Ti(x, r], A), 2(m — 2)|x|] we have 



u{t, x) — X 



ci(x) 



tP 



< \u{t, x) — voit, x) — Af i(t, x) \ + \vQ{t, x)| + A 
|ci(x)| 



Vi(t,x) 



Ci(x) 



tp 



< 3r/A- 



tP 



that is, for p = k, 



u{t,x) = A 



Cl[X) 



That gives a precise information about the time-decay of u{t, x) and shows that the estimate 
in theorem [U is optimal (for |a;|). 

In the case p = m — 1 < k the decay rate of u is determined by the decay of (long range) 
initial data and all f„ G L^p- Analogously, if we can show that vo{t,x) = co{x)t~P ^ for 
t 3> 1 then we can bound all higher perturbation orders for sufficiently small A and big t by 
the same expression multiplied by an arbitrarily small r/. To this aim we use again (|TT1) 

Cp,fcA Cm ■ (/o + /i + 5'o) 



and bound Ao(t,x) by r/ |co(x)|t^P as above. It leads to 



=:Ao(t,x) V(t,a;)G 



M(t, x) — 



< \u{t,x) - VQ{t,x) \ + 



t;o(t,a;) - 



Co (a;) 



tP 



< 2ri 



\co{x)\ 

tp 



what for p = m — 1 gives 



u{t, x) 



im-l 



That again gives a precise information about the time-decay of u{t,x) and shows that the 
estimate in theorem [1] is optimal (for t ^ \x\). 



III. NONLINEAR CASE WITHOUT THE POTENTIAL TERM 



Now, we consider a nonlinear wave equation of the form 

Du = F{u) (14) 

subject to initial data (/,(?) satisfying ([5]) with /o,/i,5'o < The nonlinear term obeys 
|F(m)| < -Fi|m|*' for |m| < 1 and \F{u) — F{y)\ < F2\u — v\ma.x{\u\,\v\)P~^. The second 
condition is satisfied e.g. for F{u) = uP with F2 = p or for F & such that |-F'(m)| < 
F2\u\P~^ for |m| < 1. 



A. Iteration 

We define an iteration scheme 

Uq := 0, 

Un+l '■= Iq 

(/,(?) + Lo(F(mO), n>0. (15) 

For it we have the following 
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Theorem 2. With f, g and F{u) as above for any m > 3, p > 1 + -\/2 and sufficiently small 
e the sequence Un converges (in norm) in Lf^^ for q = min(p — 1, m — 1) to the solution u of 
the equation [T^ . The limit u := lim„^ooW„ satisfies 

with some positive constant C depending only on p, m and e. 

Proof. For g,'Vf G L'^ and / G L'^_i with m > 3 from lemma [1] we get Ui = Io{f,g) G 
-^i?m-i- Next, if some m„ G Lf^^ with some g > 1 then, since Lq is a positive operator^, we 
have \Lo{F{un))\ < FiLQ{\un\'P) and from lemma [3] we get 

WLoiFiuML^, < Fi||Lo(|n„n|U^^ < F,C\\u4l^^, 

and hence Lo(F(u„)) G when g < p - 1. Then, Un+i G n L^^ = L^^ for 

g := min(m — l,p — 1). Hence, by induction we obtain Un G Lf^^ for every n = 0, 1, 2, ... and 

IkilU- < ||/o(/,^7)|U- < C^(/o + /i + ^7o) < 3C^£ 

< ||/o(/,^?)||l- + ||Lo(FK))|U5»^ < 3C^e + F.CWuJl^^. 
Choose e > such that FiC(6C^)%p-i < SC^ ■ min(l, 2F1/F2). Then, 



1,9 



Mr 



hence ||Mn||Lf= < ^C^^ for all > 1. As next, we show convergence of the sequence Un by 
demonstrating that it is Cauchy. 

< F2||Lo(|u„ - max(|u„|, ||l- 

< -^2^11 |m„ - M„_i| max(|M„|, |u„_i|)^"^||loo^ 



00 



with 5 := F2C{QCmsY ^ < 1, hence the iteration is a contraction in the normed space L\ 
and Un is a Cauchy sequence, because 

||?i„+l -MnlUjf^ < -Mo||ljo_^ = 5"||/o(/,5')||lj^^ < 5"3Cm£: (16) 

and for any n' > n 

n'—n—l n'—n—l 

\\Un' - UnWlf^^ < ^ \\Un+j+l - Un+j\\Lf^^ < ^ 5^^" 3Cm£: < 3(7^6. (17) 

j=0 j=0 



In fact Lq = D ^ is a measure on M^"^'^ and therefore has a positive kernel. Then, La{F) > if > 0. 
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Since is Banach, m„ has a limit u G Lf^ satisfying 

u = W.9) + Lo{F{u)) (18) 
and solving the wave equation f|T^ with the initial data ([5]). Its L^g-norm satisfies 

||m|U- < QCmS. (19) 

□ 

From ( ITTIl it follows, in the limit n' — oo an error bound 

A" 

\\u - MnlUjo^ < ^ -iC^e < Ce^^-^y''+^ (20) 

for small e. 

B. Perturbation series 

In order to be able to construct a well-defined perturbation scheme to all orders we have 
to assume additionally that F{u) is analytic at m = 0, its Taylor series starts at power p > 3 
and has convergence radius Rp > 0. Then, for small initial data 

{u,um = i6f,6g) (21) 

we introduce a perturbation series for representing the solution of f[T^ 



oo 

U = 

n=l 



Y^e^vr.. (22) 



After inserting it into (fT4|) and collecting terms according to powers of e we obtain the 
following perturbation scheme 

□^^1 = 0, {v,,v,m = {f,g) ^ v, = Io{f,g) (23) 

DVn+l = Fn{Vi,...,Vn), (^^n+l,^'n+l)(0) = (0,0) ^ = Lq (F„ (tJi , . . . , t;„) ) , (24) 

for n > 1, where Fn result from collecting the nonlinear terms with the same powers of e 

71,1 n,n 

F„(t;i,...,t;„) = 2^a>"* ■■■t;^ , (25) 

k 

1 n,m ^ ■f\j r n,m , j v-~^" n,m \ r i 

where e satisfy 2_^jn=i '^^k = n+1 and i^^^i — P every n, k. 

We call this expansion a "zero background" case because the zero-order term vq is absent. 
If a Vq term were present in the series above (i.e. the summation started at n = 0), we would 
have an additional equation Ovq = F{vq) which is truly nonlinear (opposite to the above 
system of hnear wave equations with source terms). Its solution vq represents a "background" 
around which the perturbations f„ are calculated. 

Below we show that the perturbation series converges to the solution u of the nonlinear 
wave equation (fill) and has a positive convergence radius. 
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Theorem 3. With f, g and F{u) as above for any m > 3, p > 1 + \f2 and sufficiently small 
e the series defined in ^2^-^^2^ converges (in norm) in for q = min(p — l,m — 1) to 
the solution of the equation with initial data l[21\) . 

Proof. For g,Vf G L'^ and / G L'^_^ with m > 3 from lemma [T] we get Vi = lQ{f,g) G 
L^m-i With 

< C„(/o + /i + fl-o) =: D <oo. 

Next, we prove by induction Lf'^ bounds for aU > 1 with some g > 1. Assume that for a 
given n > 1 we have Vm G for aU m < n. Then, using fl25|) . we get 



n,l n,n 

■ --Vn^ I 



Observe, that \vy 



and V \v^^ ■ ■ ■ v'^^ 



G Lf, because 



of the following estimate for 6i + . . . + 6„ = i? > 1 and 6^. > 0, m = 1, 



.n 



\w 



bi 



.bn 



<\\{t+\x\f'"'{t-\x\)^''-^^'^'^w\'\\Lo 
= ||(t+|x|)l/^(t-|x|)(^-l)/V|li^oo-- 

<\\{t+\x\){t-\x\)^'^-^^w4% 



+ \x\){t- \x\) 



\bn 

n|lL?°„' 



1,9 



which used for '■= oil'^ /p with B := oil'^ /p + . . . + al'^' /p > 1 gives 



n,l 



71. n 

I 



r oc 
-^1,9 



1^1 I 



71,1 I n.n I 



< 



"1,9 



ri.l / 

Qfc' /p 

r oc 
-^1,9 



n.n I 



< oo. 



Then, for g < p — 1 we can use lemma [3] with u : 

k 

< Ictfel • ||^i||l°° " " " ll^nllLf^ 



ri.n 



J oo 

-'"1,. 



■ ■ ■ Vn'' I to obtain 



n,l 



"1,9 



n.n 
fc 



Unfortunately, we were not able to find an estimate for |a^| being good enough to prove 
a geometric growth of ||fn|Uf=^ and guaranteeing convergence of the series (l22l) . If one tries 
so, e.g. assuming ||fm||Lf < -D™ for all m < n, then 



k ' k 



71.1 , , 7l,7l\ 

■ +...+na^' ) 



12 



The best estimate we were able to find is |a^| < Cn^ (imposing further assumptions 
on F[u)) which does not allow to close the induction argument. Therefore, we choose a 
different way and use some trick, relating the wave equation to an algebraic one. 
To this goal, we need to relate the coefficients of the power series for F{u) 



F{u) = J2bn 



n=p 



which converges for \u\ < Rp, to the expansion coefficients which result from a formal 
insertion of the series u = YlT=i^'''^k into F{u) 



\A;=1 / n=p— 1 n=p— 1 k 

By some manipulation of sums we obtain 



n,l I I n,n 

n 7 / ~r • • • "T Oif, 



where the symbol in delimiters represents the multinomial coefficient. Since there is an 
analogous relation between the absolute values of the coefficients 



n,l I I n,n 

a^' + . . . + a^' 

n.l I i„">" I n,l n,n 



\a2\ = \b-- ---1'"'= 



we observe that the series (126|) with replaced by \a^\ gives rise to a new function F 



71,1 n,n 



n=p—l k \fe=l 



such that 

oo 

F{u) = Y\bnK. (27) 



n=p 



F{u) is also analytic at m = and the convergence radius is the same as that of F{u), i.e. 
Rp = Rp what follows from standard theory of analytic functions. 
Now, instead of the system of estimates 

\\v,\\l^^<D, (28) 

n,l rL,n 

< cE Kl ■ ll^ill?- • ■ ■ IKfi^ (29) 



with q := min(m — l,p — 1), we consider a system of equations 

wi = D, (30) 

71,1 71,71 



13 



and it is easy to see (e.g. by induction) that ||fn||Lj° < w„ for all n > 1. Now comes the 
trick. Using the above relations we can find that this system is equivalent to 

CXD / OO \ 

J2 ^""wn = CF\y^ e'^Wn + De. (32) 

n=l \n=l / 

Introducing W = Yl'^=i ^"^n we can write 

W = CF{W) + De. (33) 
Since F{W) is analytic at W = 0, so is G{W) 

W - CF 

G{W) := 



D 

and also its inverse G^^{e) at e = 0, because G'(0) = 1/D > (see e.g. (real) analytic 
inverse function theorem in [5]), what follows from the fact that the Taylor series for F 
starts (as that for F) at the power at least p > 2. Then G~^{e) has a Taylor series with a 
positive convergence radius Rq-^ > 0. The solution W{e) of (!33l) can be then represented 
by a convergent series for \e\ < Rg-^ 

OO 

W{e)=G-\e) = J2^'"^n. (34) 

n=l 

In order to guarantee that this series can act as a good argument of F we choose a possibly 
smaller radius R < Rq-^ such that |W^(£)| < Rf for all \e\ < R. Then F{W{e)) can be 
represented by a convergent series (1271) in W{e). Finally, this allows us to insert this series 
into (155]) and obtain first (152]) and then the system (15U]) - (15T]) . 

Essential for the trick is that the series in (!M]) converges for all \e\ < R. Now, 
since ilVniUf^ < Wn for all n > 1, we get from the comparison criterion that the series 

X]^;^ ^^llfnlUfg converges as well for all \e\ < R. Thus, the series (122]) converges in norm 
in L?l for all \e\ < R to some u G L?l which satisfies 

OO OO 

u := J^e^Vn = eIo{f,g) + Y,^^^'L,{F^{v^, ...,v^)) 

n=l n=l 

= eh{f,g) + L, [y.^. 

= Io{ef,eg) + Lo{F{u)) 

what is equivalent to the wave equation (]T4]) with initial data (|2T]) . Uniqueness of solutions 
follows easily from theorem [21 

An important consequence of the convergence of Xl^i ^"ll'^nlUf^ is that there exist con- 
stants < M < CXD and R~^ < p < oo such that ||fn||Lf=^ < Mp" for every n > 1. □ 

Since the introduction of the auxiliary parameter e in the series expansion is only a way 
to generate the system of linear equations equivalent to the original nonlinear equation, we 
can now remove the parameter e and replace the condition on the initial data by requiring 
fo^fi^go < R- If we solve the system (|23]) - (|2^ then we obtain a solution of the nonlinear 
wave equation (fT4]) by summing up the convergent series Yl^=i '^n = u. 
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C. Optimal decay estimate 

In the nonlinear case, the iteration sequence Un is different than the perturbation sequence 
Un '■= Ylm=i'^n, therefore the question whether information about the decay rate of u can 
be read-off from the low order terms must be studied separately for both cases. On the one 
hand, in the iterative scheme the form of the source terms F{un) is much simpler than that 
in the perturbative scheme, Fn{vi, On the other hand, in practise, it is much easier 

to calculate f„'s than u^s. Below, we address both situations. 

Analogously like for the linear equation, we will have two cases depending on whether m 
is smaller or bigger than p. In the first case, the initial data will dominate the late-time decay 
rate of u, in the second case the power p of the nonlinearity, through nonlinear scattering, 
will determine the the decay rate of u. 



1. Iteration 



In analogy to the linear case, basing on a decay information for some low order term in 
the iteration sequence and on its error bound we find the exact decay rate of u. From the 
error bound (120|) it follows for large t 



\u{t,x) - Un{t,x)\ < 



{t+ \x\){t - 



< 



Cn{x)e 



{p-l)n+l 



where q := min(p — 1, m — 1). If we are able to show that some Un{t, x) = edn{x)t ^ for 
large t (the asymptotic approximation is to be understood in the following sense: 



rio>0^0<ri<rio 



3T<ooVi>r 



Un{t,x) 



edJx) 



< 1]- 



edJx) 



(35) 



i.e. the relative error t] becomes arbitrarily small for sufficiently big t, cf. ( fT2l) ). then already 
Un shows the correct decay rate, identical with this of u, because then, choosing r] := e^P~^^", 
we get 



u{t, x) — 6 



dJx) 



< \u{t, x) — Un{t, x) I -|- 



Un{t,x) - e 



dn{x) 



t1 



Cn{x) dn\X) 

< !]€ +r]e- 



for sufficiently small e. Hence the decay rate of u at late times is exactly t ^. 
In case when m>p = g + l, we have ui = /o(/, g) € L'^m-i 

\u,{t x)\ < ^"--(/o + Zi+^o) < ■ (/o + /i + ^?o) 



(t+ \x\){t- |a;|)™-2 



m—l 



for t > 2(m — 2)\x\, hence it decays faster than u and it cannot be shown that ui{t,x) = 
edi{x)t~^. It is expected that it will be true for U2 = ed2{x)t~'^, what means, that already 
U2 would have the same rate of decay as u (see jsf for such results in spherical symmetry). 

In case when m < p, we have q := m — 1. Then it should be possible to show Ui{t, x) = 
Ci{x)t~^"^~^\ what means, that already Ui would have the same rate of decay as u. 
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2. Perturbation series 



The perturbation scheme fl23l) - fl2^ can be written as 

vi = lo{f,g) 

V2=V3 = ... = Vp-i = 

Vp = Lo{Fp_i{vi, Vp^i)) = aoLo{{viY) 



Vn+l = Lo{Fn{Vi, Vn)), 



n > p. 



(36) 
(37) 
(38) 
(39) 



Assume we are in the more interesting case m > p where the tail results from the nonlinear 
scattering. Then, Vi = Io{f,g) G Lf^m-i ^^nd f„ G Lf^p-i for n > 2. If we can show that 
aoLo{{I{f, g)y) = dp{x)t~^P^^\ then already Vp has the correct decay rate, identical with 
this of u. To prove it, we need to show that eIo{f,g) and e^^^Lo{Fn{vi, ...,Vn)) for n > p 
are small relative to €^dp{x)t~^P~^\ 

Again, for vi = /o(/, (7) G Lf^_-^ the situation is obvious, e\vi\ = e\Io{f,g)(t,x)\ < 
Ci(a;)£:(l + and it is much smaller than e^dp{x)t~^P~^'> for sufficiently large t. 

From the convergence proof for the perturbation series we know that there exist constants 
M, p > such that ||fn||Lj° < Mp" for all n> 1. Hence, we can estimate the remainder 
of the perturbation series 



m=p+l 



+1 



< M ^ £™p"^ < ^ — < CeP+^ 



m=p+l 



ep 



for sufficiently small e < 1/p. It means that 



J2 ^"^™(t,x) 

m=p+l 



< 



CeP^ 



{t+\x\){t-\x\)P' 



< 



c{x)eP+^ 
(1 + t)P~^ 



for big t (and fixed x). Then 

\u(t,x) ~ ePvp(t,x)\ < \evi(t,x)\ 



m=p+l 



~ ^ tP^ 



P+i 



l-m—l 



hence with Vp = dp{x)t ^p and the relative error r] := e (cf. fl35p for the definition of "=") 

dp{x)eP 



u{t, x) 



tP- 



^ Ci{x)e ^ [c(x) + dp{x)]eP~^^ 



t 



For small e and big t such that p > ^^4^e p it follows 

<=' — clx) 



m—1 



tP- 



(x) 



u(t, x] 



dp{x)£P 



Thus Vp dominates the perturbation series for large times and small e and has the same 
decay rate as the full solution of the nonlinear wave equation u. 
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IV. NONLINEAR CASE WITH THE POTENTIAL TERM 

Finally, let's consider a nonlinear wave equation with potential 

Uu + Vu = F{u) (40) 

subject to initial data {f,g) satisfying ([5]) with fo,fi,go < ^- The nonlinear term F{u) is 
like in the previous section. 

A. Iteration 

1. Perturhative treatment of V 
As in the previous sections, we define an iteration 

Mo := 

{f,g)-XLo{Vu^) + Lo{F{un)), n>0 

We have the following 

Theorem 4. With /, g, V and F{u) as above for any m > 3, k > 2, p > 1+^/2, A < C^l and 
sufficiently small e the sequence Un converges (in norm) in for q = min(p — 1, k,m — 1) 
to the solution u of the equation (J^- The limit u := Xmin^ooUn satisfies 

{t + \x\){t — \x\)'i ^ 

with some positive constant C depending only on p, k, m, A and e. 

The proof is a combination of proofs of theorems [1] and [21 therefore we concentrate only 
on the points that differ. 

Proof. For g,'Vf G L'^ and / G L'^_i with m > 3 from lemma [T] we get Ui = /o(/, (7) G 
Next, for 6 := XCg^k < 1 there exists M > 3/(1 - 5) > and if u„ G L^^ with 
||wn||Lf < MCm£ for some n> 1 and q > 1 then from lemmas [l]l3] we get 



"1,9 



Ikn+ilU- < ||/o(/,5')IU?^, + AIIV'u^IIloc^^ + ||Lo(F(m„))||j:o 

< 3C^e + 6MCn,e + FiC^MCmfe'' < 00 

and hence Un+i G L'^^ if g := min(m — 1, k,p—l). By induction we obtain m„ G for every 
n>l. For e > such that FiC(MC^)%P"1 < min[(M(l -6)- 3)C„, M5(l - ^jc^Fi/Fs] 
we have 

hn+llUjo < 3Cm£ + SMCmS + (M(l -5)- 3)0^6 = MCmS, 
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hence llMnlUf^ < MCm£ for all n > 1. Analogously like in the previous proofs, we arrive at 

ll^in+l - UnWlf^^ < \\Lo{V{Un - Mn-l))||LJf^ + \\Lo{F{Un) - F(m„_i) ) || ^oo^ 

< XCg^kWUn - Mn-llUj", + F2C{MCmey'''^\\Un - M„_i||ioo_^ 

< 6\\Un - Un-l\\Lf;^ + 5(1 - 5)||m„ - U„,-l\\Lf;^ = S'\\Un - M„-i||l- , 

where 6' := 26 — 6"^ < 1. It follows that u„ is a Cauchy sequence (see the above proofs) in 
the Banach space and hence m„ has a limit u G L'^^ satisfying 

u = W, g) + LoiVun) + L,{F{u)) (41) 

and solving the wave equation fHOl) with the initial data ([5]). Its L^g-norm satisfies 

< MCrnS (42) 

with some (finite) constant M>3/(1— 5)>0. □ 
Moreover, by analogous considerations like in the proof of theorem [21 we find for n' > n 

1,9 ]_ — 

and in the limit n' oo 

\\u - UnWh^^^ < ^ _ (43) 

2. Non-perturhative treatment of V 

Building on the above results we can also define an alternative iteration scheme 

Uq := 

Un+i ■■= Ivif,g) + LviF{un)), n>0 

which is based on inversion of the operator □ + XV. According to the discussion in the 
introduction, it is equivalent to 

(/,(?) + Lo(F(u„))-ALo(yn„+i), n > 0. 

It converges under the same conditions as in theorem |H The proof has the only difference 
that now we have 

\\Un+l\\L^^ < ||/o(/,fi')||L- + \\Lo{F{Un))\\Lf^^ + \\\VUn+l\\Lf^^ 

< Cmifo + /l + 9o) + -^iC'll^nllLf^ + '^C'^^fc || || 

< {l-6)MCme + 6\\un+i\\Lf;^ 

what gives 

(1 - 5)MCme 
1,9 I — d 
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B. Perturbation series 



Definig a perturbation scheme for the nonhnear wave equation with potential fj40l) 

oo 

« = J]£X (44) 

n=l 

one encounters the problem of two scales which are introduced by parameters A measuring 
the strength of the potential and e measuring the strength of the initial data. There- 
fore, we propose two ways of looking at the problem: in first, we treat the potential non- 
perturbatively, in second, we assign to A a scale of some power of e. 



1. Non-perturbative treatment ofV ('A ~ e^) 

In this perturbation scheme we invert the operator □ + AV^, thus treating in a non- 
perturbative way. For the sequence Vn defined by 

V, := Iv{f,g) = W,9) - LoiVv,) (45) 

Vn+l := Lv{Fn{Vi, ...,Vn)) = Lo{Fn{Vi, ...,Vn)) - LoiVvn+i), n>l (46) 

we have the following 

Theorem 5. With /, g, V and F{u) as above for any 'm>?),k>2,p>l + \/2, A < 
and sufficiently small e the series defined in ^^^-^^D^ converges (in norm) in L'^^ for q = 
min(p — 1, k,m — 1) to the solution of the equation ( f^OP with initial data ^21\). 

Proof. The proof is essentially identical with this of theorem [3] with the following differences. 
For q := min(m — 1, k,p — 1) we obtain 

< -D + 5||t'i||L°° , , 

II II l,m — 1 'I " l.m — 1 ' 

where 6 := XCg^k and hence 



k;i roc < fl roo < — 



D 



-5 

The same modification regards all other inequalities 



< oo. 



I ^ I II II ^ II II ^ A] I 1 1 II 

Pfcl' IPi|Il°° ' ' ' \\'^n\\L°° + ^ IPn+l ||l?°„ 

l,q l,q -"-'y 

k 

which leads to 

k 

Repeating the trick used in the proof of theorem[3], we can relate this problem to the algebraic 
equation, which now becomes 

= ^ F{W) + ^^e. (47) 
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Since G(W) given by 

is again analytic, so is W{e) = G~^{e), because G'(0) = {1 — 6)/D > 0. Repeating the 
reasoning, we arrive at the conclusion that £"'||'y„||Lf^ has a positive radius of conver- 

gence. It follows that the series (jUj) converges in norm in for all e < i? to the solution 
of (HOj) with initial data (!2T|) . Uniqueness follows easily from theorem |H □ 

We can, again, remove the auxiliary parameter e and replace the condition on the initial 
data by /o,/i,5'o < R- Then, the series Yl'^=i'^n defined by (H5l) - (H6l) converges to the 
solution of the nonlinear wave equation pUj) . 



2. Perturbative treatment ofV ~ j 

If we assume that the small scale of the potential's strength A is related to the small scale 
of the initial data, say A = e"'X with a G N+, then the power series Ansatz 

oo 

M = ^£X (48) 

n=l 

inserted into the wave equation fHOj) gives 

v-n := 0, n > (49) 
Vi:=Io{f,g) (50) 

Vn+l ■= -\Lo(yVn+l^a) + Lo{Fn{Vi, .■.,Vn)), n > 1. (51) 

This system is much more appropriate for numerical techniques, because the equation on 
Vn+l is explicit, in contrast to the previous scheme, which includes implicit equations for Vn+i 
(i.e. appearing on both sides). Moreover, if we choose a := p — 1, then the lowest nontrivial 
order, Vp (all lower orders satisfy f „ = for 1 < n < p), contains both contributions from V 
and F and can be used as a good approximation to u (assuming the series converges), what 
will be discussed in the next section. 

In this case we also have a convergence result 

Theorem 6. With /, g, V and F{u) as above for any m>?),k>2,p>l + \/2, A < C~l 
and sufficiently small e the series defined in converges (in norm) in for q = 

min(p — 1, k,m — 1) to the solution of the equation [JW with initial data [21\} . 

Proof. The proof is again analogous to that of theorems [3] and with the following differ- 
ences. We have for q := min(m — 1, k,p — 1) 

ll^^lIU- < < C^mifo + fl+ go) =■ D < OO 

and 

ll^n+lIU- < \\Lo{Fn{Vi, ...,Vn))\\Lf;^ + \\\Lo{VVn+l-a)\\L^^ 

En , 1 n,7i ^ 

kfel ■ II^iIIl^ ■■■ll^nllL^ + 5||t;„+l_a|U?= • 
l,g l.g -"-'y 

k 
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with S := XCq^k- The corresponding algebraic equation becomes now 

W = C F{W) + de^W + De. (52) 
It cannot be rewritten, hke before, as G{W) = e, but it can be written as 

G{W, e):=W ~C F{W) - 6e^W -De = 0. 

Since a is a positive integer number, G{W, e) is analytic in both variables around the point 

G(0,0) = 0. Moreover, dG{W,0) /dW\w=o = l-Se" = 1- XC^^k > 0. Then, by (real) 
analytic implicit function thorem (see e.g. [5|]), there exists a unique function W{e) such 
that G{W{e),e) = 0. Then W{e) has a Taylor series representation with positive radius of 
convergence. Repeating the reasoning of the previous proofs, we arrive at the conclusion 
that ^"ll'^nlUfq ^ positive radius of convergence R> 0. It follows that the series 

converges in norm in for all < e < R to the solution of (140 p with initial data 
(pTj) . Uniqueness follows again from theorem |H □ 



C. Optimal decay estimate 

1. Iteration with perturbative treatment ofV 

From ( l43l) we have 

5"(2-5)" 3Cm€ c„(x)A"e 

\u{t,x) - Un{t,x)\ < -Jj^ZTsy- ■ (t+|a;|)(t-|x|)^-i " (l + t)" 

for sufficiently small e and A (such that XCq^k = S < 6o < 1). If we are able to show 
that some Unit.x) = edn{x)t~'^ ^ for large t (the asymptotic approximation "=" is to be 
understood in the sense defined in fl35|) . with the relative error rf) then already Un shows the 
correct decay rate, identical with this of m, because then, choosing rj := A", we get 



u{t, X) — 



t1 



^ \dn{x) + c„(3:)]A"g 
t<i 



For small A it follows 

Again, in case when m > p, we have Ui = /o(/, g) G L'^m-i 

I Cm-{fo + fi + go) ^2a^-(/o + A + 



(t+ |a;|)(t- - (l + t)"^-i 

for t > 2(m — 2)\x\, hence it decays faster than u and it cannot be shown that Ui{t,x) = 
edi{x)t~'^. It is expected that it will be true for U2 = ed2{x)t~'^, what means, that already 
U2 would have the same rate of decay as u (see ^ for such results in spherical symmetry). 

In case when m < p, we have q := m — 1 . Then it should be possible to show ui{t, x) = 
di(x)t-('"-i), what means, that already ui would have the same rate of decay as u. 
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2. Perturbation series with perturbative treatment of V 

Consider the system fH9l) - fl5Tl) and choose the constant a := p — 1 so that V2, Wp-i = 
and at the order Vp both effects, the nonhnear and hnear (potential) scattering, appear 
simultaneously 

v-n := 0, n > (53) 
vi = Ioif,g) (54) 

V2 = V3 = ... = Vp_i = (55) 

Vp = -XLoiVvi) + LoiFp^iivi, ...,Vp.i)) = -XLoiVvi) + a^L^{{viY) (56) 

Vn+l = -XLo(yVn-p+2) + Lo{Fn{Vi, f„)), U > p. (57) 

Consider only the more interesting case m — 1 > min(p — 1, k) =: q. If we can show that 
Vp = dp{x)t~'^, then already Vp has the correct decay rate, identical with this of u. To prove 
it, we can repeat the reasoning from the section when we treated nonlinear wave equation 
without the potential term, because the only fact, which we use is that the perturbation 
series "^^=1 ^""^n has a positive radius of convergence and this is here guaranteed by theorem 
[6l Analogously, we obtain 

^ rfp(x)£P 

for alH > T and sufficiently big T = T{6), so Vp dominates the perturbation series for large 
times and small e and has the same decay rate as the full solution of the nonlinear wave 
equation u. 

This is the simplest setting for applications. Here, we only need to solve (approximately) 
two linear wave equations, ( !54|) and (l56| . in order to determine the decay rate for solutions 
of ( HOl) . This is the starting point of [3| where we solve the two equations under spherical 
symmetry. 



APPENDIX A: SOME USEFUL ESTIMATES 

The first two lemmas we cite from 0]. 
Lemma 1. Let the data {f,g) G L'^_i x L'^ with m > 3 satisfy 

fo ■= ll/IUjS_i < /i := I|V/||l- < 9o ■= IklU- < oo. 

Then there exists a unique weak solution v{t,x) = /o(/, (?) of the free wave equation 

Uv = 0, v{0, x) = f{x), dtv{0, x) = g{x) 

which satisfies 

\\v\\l^^_, < C{f,g) := ■ {go + fi + fo)- 
Lemma 2. Let the source F satisfy for some q > 2 and 1 < p < q 

Fo:= < oo. 
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Then there exists a weak solution v{t,x) = Lq{F) of the free wave equation with source 

Dv = F, 

and null initial data v{0,x) = 0, dtv{0,x) = 0. Moreover, it satisfies 

< Cp^qFo. 

Next lemma we cite after Asakura [1, Cor. 2.4 and Eq. 2.33] and state in our notation. 
Lemma 3. Let u G C^(M+"*'^) fl for some q> I. Then for any p > I + \/2 

with some C > provided q < p — 1. 

Note, that it is a consequence of lemmaO but only when p > 3, while for 1 + ^/2 < p < 3 it 
requires a more general proof. It can be easily deduced, also for weak solutions u E C°(M+^^), 
from a more general estimate 



Lemma 4. // 

\F{t,x)\< ^ 



{t+ \x\)P{t - \x\)i 



with p > 2, q > 1 then 

\LoiF)it,x)\< ^ 



{t+ \x\){t- 

with some positive constant C . 
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